Optical absorption from solvation-induced polarons on nanotubes 
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When an excess charge carrier is added to a one-dimensional (ID) semiconductor immersed in a 
polar solvent, the carrier can undergo self-localization into a large-radius adiabatic polaron. We ex- 
plore the local optical absorption from the ground state of ID polarons using a simplified theoretical 
model for small-diameter tubular structures. It is found that about 90% of the absorption strength 
is contained in the transition to the second lowest-energy localized electronic level formed in the 
polarization potential well, with the equilibrium transition energy larger than the binding energy of 
the polaron. Thermal fluctuations, however, cause a very substantial - an order of magnitude larger 
than the thermal energy - broadening of the transition. The resulting broad absorption feature may 
serve as a signature for the optical detection of solvated charge carriers. 

PACS numbers: 78.67.-n, 71.38.-k, 31.70.Dk 



I. INTRODUCTION 

One-dimensional (ID) semiconductor (SC) nanostruc- 
tures (e. g., nanotubes and nanowires) in contact with 
polar solvents is an interesting class of systems of par- 
ticular relevance to applications and processes involving 
fundamental redox reactions [H, [3i When an excess 
charge carrier (an electron or a hole) is added to such 
a structure, the carrier can become solvated resulting in 
a polaron, a self-consistent combination of a localized 
electronic state and a dielectric polarization pattern of 
the surrounding medium d, H, d, 0, [1| ■ The long-range 
Coulomb mechanism of the polaron formation here is 
analogous to the well-known three-dimensional (3D) po- 
larons in polar SCs |9|. [Tol . [T]| and solvated electrons in 
polar liquids [Tll.[T^.ll3j|. Distinct from those cases, how- 
ever, are the confinement of the electron motion in a ID 
nanostructure and its structural separation from the 3D 
polarizable medium. 

The physical properties of polarons are quite differ- 
ent from the band electrons frequently discussed in the 
context of electronic transport in nanotubes [ij, and 
nanowires p^ . We have already emphasized the ener- 
getic significance of the solvation by finding [^,@] that the 
binding energy of the resulting ID polarons could reach a 
substantial fraction, roughly one third, of the binding en- 
ergy of Wannier-Mott excitons, the well-known primary 
photoexcitations in many ID SCs. This may lead to en- 
hanced charge separation. On the other hand, the mobil- 
ity of solvated charge carriers is drastically reduced due 
to the dissipative drag of the medium [1, |^ [l3| • In this 
paper we address qualitative theoretical expectations for 
the local optical absorption from solvation-induced ID 
polarons. This kind of the optical absorption has been a 
powerful tool in optical detection of electron-lattice po- 
larons in conjugated polymers as well as of 3D solvated 
electrons. 

While the polaronic effect and the features we dis- 
cuss have a generic character, specific illustrative calcu- 
lations in this paper are done using a simplified model 



representation for small-diameter tubular structures. A 
widely known example of such structures is semicon- 
ducting single- wall carbon nanotubes (SWCNTs). We 
note that redox chemistry of CNTs has been deemed an 
"emerging field of nanoscience" [l^ and solvatochromic 
effects in CNTs are being intensely researched . 

As discussed in more detail later. Fig. [T] displays cal- 
culated functional dependences for the binding energy, 
E\j, and the spatial extent, Z, of ID adiabatic polarons in 
terms of the tube radius R and convenient scales of the 
Bohr radius and Rydberg energy for the corresponding 
3D Coulomb problem [36| : 



as 



e*fi^/mq^, Ky — /2e*aB- 



(1) 



Here q is the carrier charge, m its effective mass, and 
e* the effective dielectric constant of the uniform 3D 
medium. The length I has been defined as the fuU-width- 
at-half-maximum (FWHM) of the localized ID electron 
charge density. In what follows we refer to the ratio 
as/R as the confinement parameter. It is well known 
that the spatial confinement generally results in amplifi- 
cation of the Coulomb binding ^20*1, and Fig. [T] illustrates 
both the magnitude and the growth of Eb/Ky for ID po- 
larons with the increased confinement parameter. For a 
proper perspective, one should compare ^ those results 
with the polaron binding energy in higher dimensions, 
where the classic Pekar's result [21|, i22| for Ej, in 3D is 
only ~ 0.1 Ry, while in 2D systems [2j| it increases to 
~ 0.4 Ry. Our variational calculations @ have shown 
that the transition to a ID polaron structure (charge 
uniformly distributed over the tube circumference) takes 
place at as / i? > 1 . 

The effective dielectric constant in Eq. ([T]) varies for 
different solvents and can be determined from the well- 
known d, [nl, m, [H relationship 

1/e* = I/eco - l/e^; 

typical [l^l static and high-frequency Cqo constants 
satisfy eg 3> Coo so that e* ~ Eoo. For SWCNTs, effec- 
tive band masses m of carriers generally depend both on 
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FIG. 1: The binding energy, Eb, and the linear size, I, of ID 
adiabatic polarons as a function of the confinement parameter 
as/R- The size I is defined here as the FWHM of the local- 
ized electron density. Different curves correspond to different 
dielectric conditions, as explained more fully in the text: filled 
(solid lines), hollow (long-dash), and polarizable (short-dash) 
nanotube screening models. 



tube radius and chirality; see, e.g., Ref. [TH for a com- 
pilation of some theoretical data. If, for instance, one 
were to choose a set of representative numerical parame- 
ters: m = 0.05TOe {me being the free electron mass) and 
e* = 3, then Eq. Q results in as = 32 A, Ry = 76 meV, 
and, with a roughly estimated as/R ~ 4, ~ 0.1 eV. 
Along with results demonstrating overall scaling with the 
confinement parameter, we will be using specific values 
of as/R — 4 and 10 for some plots. One should under- 
stand that those are meant to serve mainly illustrative 
purposes. 

Figure [5] displays an example of the equilibrium ID 
polarization potential well calculated for as/i? = 4. 
The corresponding pattern of the medium polarization 
is self-consistently stabilized by the electric field of a 
charge carrier occupying the lowest-energy localized elec- 
tronic level formed in this well (the ground state of the 
polaron). Along with the lowest, also shown are the 
two next higher-energy empty levels. The results dis- 
cussed below indicate that it is the transition between 
the two lowest-energy levels that has the largest oscil- 
lator strength and dominates the optical absorption in 
the ground state. The long-range behavior of the poten- 
tial well is Coulombic, and hence it contains many more 
higher-energy localized states (not shown in Fig. ^ con- 
verging to the onset of the continuum spectrum at zero 
energy in the figure. 

Our consideration in this paper assumes that the po- 
laron binding energy Ei, is sufficiently larger than the 
thermal energy fc^T so that the thermal occupation of 
higher-energy electronic states can be neglected. Even 



FIG. 2: The equilibrium self-consistent ID potential well (a; is 
a coordinate along the tube axis) along with the three lowest- 
energy electronic energy levels calculated for ub / R ~ 4:. The 
optical transition with the largest oscillator strength is indi- 
cated by the arrow, the dashed lines showing the localized 
electronic wave functions of the participating states. 

so, thermal fluctuations are found to very significantly 
broaden the optical transition i n q uestion, as is also the 
case for 3D solvated electrons [23]. In the presence of 
the polaron, the fluctuations involve combined dynamics 
of the electronic charge density and the solvent dielec- 
tric polarization. A major contributor to the broadening 
is the resulting ID polaron "breathing" mode, qualita- 
tively corresponding to fluctuations of the polaron length. 
Without pursuing the exact calculation of the absorption 
line shape, we will provide an approximate estimation of 
the broadening due to the breathing mode within the 
classical Franck-Condon framework 26]. The described 
broad absorption feature could be used for detection of 
the solvation-induced ID polarons and an estimate of 
their binding energy in prospective experimental studies. 

II. THE GROUND-STATE AND OPTICAL 
TRANSITIONS IN EQUILIBRIUM 

In the single-particle ID continuum adiabatic frame- 
work, the excess charge carrier is described by the wave 
function ipix); it responds to the electrostatic potential 
(j){x) via the corresponding Schrodinger equation: 

- + <iH^)M^) = E^M^)- (2) 

In the ground state, the charge carrier occupies the lowest 
energy {i = 0) state so that the resulting ID electronic 
charge density is 

p{x)^q\M^)\^- (3) 
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At equilibrium, the self-consistency requires that this 
charge density stabilizes the very electrostatic potential 
used in Eq. 



^(a;) = / dx'G{x - x')pix'), 



(4) 



where kernel G is the appropriate electrostatic response 
function assuming the translational invariance along the 
X-axis. Solving Eqs. ©-([I]) together yields the equilib- 
rium ID polaron distributions that we will be denoting as 
(/)(a;) and p{x). One notes that the equations determining 
the equilibrium correspond to the minimum of the total 
adiabatic energy of the system: 



77'tOt TP 

^0 — -^0 



Ep, 



(5) 



consisting both of the electronic energy Eq in Eq. ^ 
and the energy stored in the dielectric polarization of the 
medium: 



Ev 



1 



dx dx' (t){x)G-\x-x')4>{x'), 



where G~^ is the kernel inverse to G (see also Ref. Q). As 
the continuum onset in Eq. ^ is set at zero energy, the 
"depth" of this minimum is the polaron binding energy: 



mm{El°'} = -Eb. 



(6) 



The data for Eb and equilibrium p{x) and q4>{x) have 
been used in Figs. [T] and [D Other adiabatic energy sur- 
faces El°*- — Ei + Ep are defined similarly to Eq. ([5]) . 

The electrostatic potential (j){x) in Eq. ^ is due to 
only the slow (orientational) polarization of the sur- 
rounding medium and hence should not include the "in- 
stantaneous" self-interaction of the charge carrier with 
itself. A proper separation of the slow polarization 
resp onse is a standard step in the polaron problem 
[m, [Hill [111- The kernel G in Eq. Q can therefore be 
written as 



G{x)=Gsix)-Gooix), 



(7) 



where indices s and oo denote that the standard elec- 
trostatic potential problem ^ is solved with dielectric 
conditions corresponding respectively to full (static, s) 
screening or to screening only by "fast" (oo) components 
of polarization. In the polaron context, "fast" means op- 
erative on time scales shorter than fi/ E^. In the model we 
consider, the continuum solvent is characterized by the 
corresponding dielectric constants and Coo- Another 
"fast" contributor to the screening is the polarizability 
of the ID SC nanostructure itself related to interband 
electronic transitions. 

In Fourier space, g{k) — J dx e~^''^ G{x), Eq. ([7|) can 
be conveniently re-written as 



= 9o{k) 



1 



1 



.(fc) 



(8) 



where the fc-dependence of the corresponding dielectric 
functions reflects a possible spatial dispersion of the 
screening due to the geometry of the system. 

The bare, unscreened, response for the ID tubular ge- 
ometry (the "elementary" charge distributions are uni- 
form rings of radius R) is given by [27| 

go{k) = 2Io{kR)Ko{kR), 

where I{x) and K{x) are the modified Bessel functions 
appearing in electrostatic problems with cylindrical sym- 
metry [28|. The general form of the corresponding k- 
dependent dielectric functions can be represented fsl] as 

e{k) = kR[ei Ii{kR)KQ{kR) + £2 IoikR)Ki{kR)]+n{k), 

where ei is the dielectric constant of the dielectric 
medium inside the tube, €2 the dielectric constant of the 
unbounded medium outside the tube, and n(fc) the con- 
tribution arising from the polarizability of the tubular 
surface itself. 

If the latter contribution is neglected and the media 
inside and outside are the same (we will be calling this 
model case "filled"), there would be no spatial dispersion 
in the screening (uniform medium) and response ([8]) takes 
a simple form of 



m^-9oik)/e* 



(9) 



In order to evaluate the magnitude of the effects arising 
from different dielectric screening conditions, we also ex- 
amine the model cases called "hollow" , where the interior 
of the tube does not contain any medium (ei = 1), and 
"polarizable" , in which the "hollow" model is augmented 
by the tube polarizability. For the model of the latter, 
we will use the following approximate representation of 
the polarizability of semiconducting SWCNTs calculated 
in Ref. fH: 



Uik) 



&.2{kRfh{kR)Ko{kR) 
l + l.Q{kRY-^ ■ 



Figure [T] compares the results for the ground state of 
the polaron calculated with the three dielectric screening 
models, in which we used representative values of £00 = 3 
and Cs = 40 for the solvent medium. It is evident from the 
comparison that in the parameter range we study, differ- 
ent screening models do not lead to drastically different 
results. The variations observed are quite understand- 
able and are largely caused by the effective changes in 
the magnitude of the fast screening, which is decreased 
by going from the filled to hollow model and is increased 
in the polarizable model. Relatively small changes are 
also observed in the results for individual electronic en- 
ergy levels displayed in Fig. [31 Panel (a) of that figure 
just shows the scaling of the three lowest electronic en- 
ergy levels Ei (i = 0, 1, 2) from Eq. ([2]) as formed in the 
equilibrium self-consistent potential well (see Fig. ^ . 

For the optical absorption from the ground state, a 
Franck-Condon framework discussion requires electronic 
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FIG. 3: (a) The energies Ei of the three lowest electronic en- 
ergy levels in the ground state equilibrium as a function of the 
confinement parameter for filled (solid), hollow (long-dash), 
and polarizable (short-dash) screening models, (b) The ratio 
of the polaron binding energy, Et, to the electronic transition 
energy, Eio, in equilibrium calculated with the same models. 



transition energies 

Em = E, - Eo. (10) 

The dominant dipole electronic transition turns out to 
be i = * = 1, and Fig. [3]Jb) shows the relationship 
between the polaron binding energy Eh and the equi- 
librium value of EiQ. (To compare, the ratio of these 
energies for 3D Pekar's polaron is about 0.8 More 
detailed information on this transition for our ID polaron 
is displayed in Fig.[4l the most important conclusion from 
which is that the transition contains nearly 90% of the 
total absorption strength across a wide range of the con- 
finement parameter. This overall dominance is similar to 
the one taking place for 3D solvated electrons [25j . 

Since Eq. ^ has a form of a standard continuum 
Schrodinger equation, the overall absorption satisfies the 
optical sum rule familiar from atomic systems [20j : 

i 

where the oscillator strength of the transition to the ith 
state 

f^O^^^\{M^)\^/^x\M^))\\ (11) 

In terms of the frequency w-dependent oscillator strength 
density [lO] , the sum rule per one polaron reads 



FIG. 4: Oscillator strength and transition energy for the op- 
tical transition between i — and i = 1 electronic eigenstates 
in the equlibrium potential well. Calculations done with dif- 
ferent screening models are denoted as in Figs. [1] and |3] 

where, for sharp discrete transitions pop . 

fi'^) = X! '^(^ ~ '^^o)' " Eio/h. (13) 

i 

We will discuss in the next section how the thermal 
broadening of these transitions occurs resulting in a con- 
tinuous f{L0). 

As Eq. pT|) clearly indicates, the optical absorption 
in our case is due to the electric field polarized parallel 
to the tube axis. The data in Fig. 2] has been derived 
by calculating matrix elements in Eq. (jlip between the 
states in the equilibrium potential well (and the sum rule 
verified in all cases). 



III. THERMAL BROADENING OF THE 
OPTICAL ABSORPTION 

In the Franck-Condon framework, the optical absorp- 
tion from the polaron is realized by virtue of "verti- 
cal" electronic transitions between the adiabatic poten- 
tial surfaces meaning that the polarization state 
of the environment does not change while the transition 
takes place. This picture of slow modulating dynamics of 
the environment leads to what is called inhomogeneous 
broadening of the absorption lineshape (see Chapter 5 of 
Ref. [26] for an extensive discussion) and is well justified 
for our application to the case of the orientational solvent 
polarization. Indeed, the parameter range of our interest 
and practical significance corresponds to the regime of 
Eb > ksT > Ti/tl, and, hence. 



\ f{L0)dL0^l, 




(12) 



TL :^h{kBTEby^/\ 
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which is a condition of the appHcabihty of the semi- 
classical Franck-Condon picture [26[. Here tl is the lon- 
gitudinal relaxation time of the solvent defining the rel- 
evant time scale for the solvent orientational dynamics 
and equal to [ll,[2l|3l| 

TL = (eoo/es)TD 

for the Debye solvent described by the w-dependent di- 
electric function 



e(c.) 



1 — iujTo 



(14) 



Typical solvents are characterized by a wide range of 
ranging from fractions to tens of ps [1, [l^l • 

Thermal fluctuations result in variations of the poten- 
tial distribution (j){x) in Eq. Q from the equilibrium pat- 
tern (^{x) discussed in the previous section, thereby mod- 
ulating electronic wave functions and transition energies 
([To)) . Averaging over various fluctuations would thus lead 
to broadening of the discrete transitions in Eq. In 
the discussed semi-classical picture, the weighting of var- 
ious fluctuations in the ground state is described by the 
Boltzmann factor associated with fluctuations of the to- 
tal adiabatic energy ([5]) from the minimum equilibrium 
value (O. Of course, fluctuations in the presence of the 
polaron are not the fluctuations of the solvent alone but 
rather combined fluctuations of the medium and localized 
electronic charge density. 

One can get an insight into the nature of thermal fluc- 
tuations of the polaron shape by using explicit results 
derived in Ref. [8] for small- amplitude fluctuations in a 
simplified model of the charge carrier interacting with the 
Debye solvent (fTi|) with response function Within 
that model, shape fluctuations arc governed by ampli- 
tudes a„(t) in the expansion of the time-dependent elec- 
tric potential 4>(x, t) and self- localized carrier charge den- 
sity p{x,t) over the normal dielectric relaxation modes: 



(j){x,t) = 0(x-e(t)) + ^a„(i)</)„(a:-C(i)), (15a) 

n>l 

p{x,t) = p(a:-e(0) + 5I««W/^«(^~^W)' (15b) 

ri>l 

where the first terms in the r. h. s. of Eqs. (|15p cor- 
respond to the equilibrium static patterns. Both equi- 
librium patterns and normal modes are centered around 
arbitrary polaron centroid positions S,(t). The dynamics 
of ^ (t) results in the diffusion of the polaron as a whole 
and does not affect the optical absorption lineshape. The 
spatial patterns of the normal modes, ((/)„(x), p„(x)), fol- 
low from the solutions of the generalized eigenvalue (A„) 
problem: 



Pn{x) = - dx'L{x,x')(f>n{x'), 



(16a) 



4>n{x) = --^ / dx'Go{x - x') pn{x'), (16b) 



where spatial kernel L is determined by the response of 
the electronic subsystem (see Ref. |^ for more detail). 
Equations (fT6|) yield shape- modulating modes (n > 1) 
with growing eigenvalues A„ > 1 as well as the zero- 
frequency translational mode with Aq = 1, which is ex- 
cluded in Eq. (fT5|) in favor of the collective coordinate 

m- 

For our purposes here, it is convenient to choose the 
normalization of the modes as 



dx pm{x) (j)n{x) = 



(17) 



note that coefficients a„(t) are then measured in units of 
the square root of energy. The equations of motion for 
so normalized modes read m 



where mode relaxation times t„ are given by 

= 1 - 1/A„ 

and perturbations ry„ by 



VniS.^t) ^ ~\n / dxpn{x)F{x -y- £,,t). 



(18) 



(19) 



(20) 



As the eigenvalues A„ grow with the mode index n, the 
relaxation times t„ in (|19p converge to the solvent longi- 
tudinal relaxation time tl (see [8| for an illustration of 
this convergence). 

For thermal dielectric fluctuations, F{x^ t) in (|20p is a 
random process with a vanishing average and correlations 
established by "bare" solvent fluctuations Q: 



{F{x,t)F{x' ,t')) 



e^TL 



Go{x-x')d{t-t'). (21) 



It then follows from and ^ that different 

modes n and m exhibit no dynamic correlations, and the 
only relevant (for the same centroid coordinate) non-zero 
correlation is 

{Vni^,t)7Jnitt'))=2DSit-t'), D = kBT/TL, 

with the effective diffusion coefficient D independent of 
the mode index with the chosen normalization. 

The random dynamics described by Eq. p8)) is thus the 
Ornstein-Uhlenbeck process [s^l, which leads, at times 
t ^ r„, to the Gaussian distribution of a„ with the vari- 
ance established by 



the probability 



P(a„) oc exp ( --jr-rj^ ) i Un{an) 



2 V A„ , 



(22) 
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FIG. 5; Modulation of the two lowest electronic energy levels 
Eq and El by the normal fluctuations modes, Eqs. (|15p . as 
a function of mode amplitudes in arbitrary units. Solid lines 
show the modulation by the even mode n = 1 and long-dash 
lines by the odd mode n = 2. The confinement parameter 
used for this illustration is as/ii = 10. 



Equation features the classic Boltzmann distribution 
for independent fluctuation modes with the effective po- 
tential energies t/„(a„). The behavior of J7„ is consistent 
with the behavior (|19p of relaxation times. 

As the normal mode index n grows, the interaction 
of the solvent polarization with the electronic density in 
the ground state generally decreases, likewise the mod- 
ulation effect on the transition energy Eiq generally di- 
minishes. The largest modulation effect occurs due to 
the even normal fluctuation mode n ~l, corresponding 
to the shallowest potential energy Ui . This modulation is 
illustrated in Fig. [S] comparing the effects of two normal 
modes on the electronic energies Eq and Ei . It is evident 
that the major contribution comes from the modulation 
of the ground state level Eq by the n = 1 mode. As dis- 
cussed in Ref. Q, the spatial pattern of this mode cor- 
responds qualitatively to the modulation of the polaron 
length Z, we will accordingly call this the "breathing" 
mode. 

While illuminating the character of the fluctuations 
and their modulation effects, the small-amplitude de- 
scription turns out to be insufficient for realistic T on 
the order of room temperature. In this paper we do not 
pursue a comprehensive analysis of the absorption line- 
shape that would then be required. Instead, we will pro- 
vide a qualitative assessment of the effect of the polaron 
breathing using a simplistic single-parameter scaling de- 
scription similar in the spirit to the approach used in a 
recent study [ssj of a continuum adiabatic model of 3D 
electron solvation. In this description, the fluctuations 
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FIG. 6: Variations of the electronic energy SEq, polarization 
energy 5Ep and the total adiabatic energy 5Eq'^ with the 
scaling parameter c. The confinement parameter used for 
this illustration is as/R = 10. 

are restricted to modulate the electrostatic potential via 
a scaling ansatz of the equilibrium pattern: 

4){x) = c4){cx) , (23) 

where c > is the scaling parameter. Its deviations from 
the value of c = 1 results in variations of all quantities 
relevant for the optical absorption. It should be noted 
that the results of small deviations from c = 1 turn out 
to be quite close to the results we derived with exact 
small- amplitude breathing fluctuations. 

Figure [5] displays variations of different energies 
5E{c) — E{c) — E{c — 1) with the scaling parameter c. 
It is evident that the shallow response of the total energy 
6Eq°^ to fluctuations around the equilibrium is a result 
of a near compensation of much larger variations of the 
polarization energy SEp and the purely electronic energy 
SEq, as indeed should be the case for the ground state 
adiabatic surface. It is this strong variability of -Bo in 
response to the polaron breathing that has already been 
noted in Fig. O and is largely responsible for the resulting 
enhanced broadening of the lineshape. The "mechanism" 
of the enhancement is also illustrated in Fig. [7] display- 
ing two adiabatic energy surfaces, Eq°^ and , between 
which the dominant optical transitions take place. The 
behavior of E\°^ here reflects large variations SEp of the 
polarization energy without a comparable compensation 
of the electronic energy that occurs for Eq°^ . Figure [7] 
shows how a population on i?g°* within a relatively small 
energetic broadening translates into a much broader 
distribution of the transition energies. 

To model the absorption lineshape as a function of 
frequency w, we work with an ensemble average over var- 
ious fluctuations represented by different values of the 
scaling parameter c in Eq. (|23p . The resulting oscillator 
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FIG. 7: Total adiabatic energies and El°^ as functions 
of the scaling parameter c calculated here for the confinement 
parameter ub/R = 10. The plot illustrates how fluctuations 
of small energy At translate into much larger variations A^j 
of the electronic transition energies. 

strength density in Eq. is then a Boltzmann average 
of Eq. taking the form of 

f{co) = / dcP{c) ho{c) 5{uj - w,o(c)), (24) 
where the probability 

-'"'--(-^)//''"-(-^) 

determines the weighting of fluctuations on the ground 
state adiabatic surface i?Q°*. Both transition frequencies 
CL'io(c) and oscillator strengths /io(c) in Eq. (I24|) are func- 
tions of the scaling parameter c. 

Figure [8] shows examples of the absorption lineshapes 
calculated with Eq. (p4|l for two values of the confine- 
ment parameter and two different values of temperature 
T . The results are broad peaks around the respective 
equilibrium transitions frequencies Eiq. The broaden- 
ing of the absorption fines, as we discussed, is strongly 
enhanced, roughly an order of magnitude larger than 
the thermal energies in our examples. In this respect, 
it is similar to the results known for 3D solvated elec- 
trons ,25]. It is this large and confinement-parameter- 
dependent broadening that is the emphasis of Fig. [H] 
rather than details of the model lineshapes. 

IV. CONCLUSIONS 

As a result of solvation by the surrounding sluggish po- 
lar medium, excess charge carriers on ID semiconductor 
nanostructures can undergo self-localization into ID po- 
larons whose properties are quite different from band car- 
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FIG. 8: The model absorption lineshapes calculated as de- 
scribed in the text for as/7? = 4 and aB/R = 10 and for two 
different temperatures: ksT — 0.2 Ry (long-dash curves) and 
ksT — 0.3 Ry (solid curves). Vertical short-dash lines show 
the positions of the respective electronic transition energies 
Eio in the equilibrium configurations. The low-frequency tail 
of f{ijj) for the case of as / R = 4 is not shown completely as 
the computational scheme used becomes inadequate. 

riers. In an effort to establish experimentally testable sig- 
natures of the polaron formation, in this paper we exam- 
ined the local optical absorption from solvation-induced 
ID polarons within a simplified single-particle adiabatic 
framework for tubular structures. 

The calculations indicate that a new broad absorption 
feature is expected, the peak energy of which is com- 
parable with (but larger than) the binding energy of the 
polaron and whose width is much larger than the thermal 
energy. This feature arises mostly due to the electronic 
transition between two lowest-energy localized electronic 
levels formed in the polarization potential well. We iden- 
tified the polaron length fluctuations ( "breathing" mode) 
as the major source of the broadening. 

Similarly to the enhanced binding of ID excitons [20j . 
the binding of ID polarons is stronger than in higher di- 
mensions, and we explored polaron structure and optical 
transitions as a function of the conflnement parameter 
a B /R, where as is the 3D Bohr radius ^ and R the ra- 
dius of the tubular structure. Of particular relevance for 
the optical absorption is the equilibrium transition en- 
ergy Eio, which approximately determines the position 
of the absorption peak (Fig. [8]). Within the range shown 
in Fig. 31 the scaling of this energy may be approximated 
as 

Sio/Ry « iae/Rf-^. (25) 

Note that the relation between Eio and the binding en- 
ergy Eb also changes with the confinement parameter 
(Fig. [3]). Together with the parametric dependence ([T]) 
of the Rydberg energy Ry, Eq. allows to assess how 
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the position of the absorption peak scales with system 
parameters (compare to scahng of excitons in SWCNTs 
pl.[35l|). It follows then that the effective dielectric con- 
stant e* of the solvent affects it as 

If the effective mass m of the carrier was independent of 
the tube radius R, then ([25]) would result in 

£^10 cx m°-3. 

In SWCNTs, however, m does depend on tube radius, 
and, for not very small radius tubes, this dependence is 
roughly m (x 1/R [H, [s^l- In this regime, the confine- 
ment parameter becomes practically independent of R, 
and the transition energy scales as 

Eio oc m oc R~^. 



Numerical estimates provided in the Introduction suggest 
that the absorption peak may be expected at energies on 
the order of 0.1 eV, subject, of course, to variations of 
the effective system parameters. 



To reiterate, illustrative results in this paper have been 
derived using simplified model considerations. More ac- 
curate treatments, particularly of temperature effects, 
may be needed for detailed comparisons with experimen- 
tal data when it becomes available. 
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nology Transfer Program. 
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